Bond-impurity induced bound states in disordered spin-1/2 ladders 
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We discuss the effect of weak bond-disorder in two-leg spin ladders on the dispersion relation of 
the elementary triplet excitations with a particular focus on the appearance of bound states in the 
spin gap. Both the cases of modified exchange couplings on the rungs and the legs of the ladder are 
analyzed. Based on a projection on the single-triplet subspace, the single-impurity and small cluster 
problems are treated analytically in the strong-coupling limit. Numerically, we study the problem 
of a single impurity in a spin ladder by exact diagonalization to obtain the low lying excitations. At 
finite concentrations and to leading order in the inter-rung coupling, we compare the spectra obtained 
from numerical diagonalization of large systems within the single-triplet subspace with the results 
of diagrammatic techniques, namely low-concentration and coherent-potential approximations. The 
contribution of small impurity clusters to the density of states is also discussed. 



I. INTRODUCTION 

Since the synthesis of spin ladder materials such as 
SrCuaOg (Ref. H or (Sr,Ca,La)i4Cu2404i (Rcf. 0) and 
the later discovery of superconductivity under high pres- 
sure in Sro.4Ca13.6Cu24O41.84 (Ref. 0), spin-1/2 two-leg 
spin ladders have been on the focus of theoretical activ- 
ities (see Refs. and |^ for review). This model is of 
particular interest, as the ground-state is nonmagnetic 
and it is an example for a spin liquid. Many properties 
of the pure system, such as the dispersion of elementary 
excitation9Sii2iS and the thermodynamicsiS are well un- 
derstood. The elementary excitations are propagating, 
massive triplet-modes. 

A natural extension of the pure spin model comprises 
the inclusion of impurities. In this paper, we discuss 
the effect of impurities on the dispersion relation of the 
elementary triplet excitations and focus on the appear- 
ance of bound states in the spin gap. Such states could 
be visible in resonant experiments. One may distin- 
guish between different kinds of impurities. First, mag- 
netic ions, such as the Cu^+ ions in SrCu203, can be 
replaced by nonmagnetic ones such as Zn (see Ref. ITTI 
and references therein), effectively removing a spin-1/2 
moment, or by other ions with the same or a different 
effective moment. The replacement of the spin-carrying 
ion will be referred to as a site impurity. Note that it 
is also conceivable that a site impurity leads to modifi- 
cations of the exchange couplings to neighboring sites. 
Second, and this is what we mainly have in mind in 
this study, the exchange paths themselves can be mod- 
ified by doping the bridging X-ions in, e.g., Cu-X-Cu 
bonds, realizing what we call a bond impurity (or sim- 
ply impurity) in the following. Such a situation is de- 
scribed in Ref. [13 for the alternating spin chain system 



(CH3)2CHNH3Cu(CUBri_^)3 where CI and Br ions are 
substituted with each other. Furthermore, a spin lad- 
der material exists, namely (C5Hi2N)2CuBr4 (Ref. IT^ 
where one could think of analogous doping experiments. 
This material is suggested to contain two-leg spin lad- 
ders in the strong-coupling limit, i.e., the coupling con- 
stant along the legs Jl is small compared to the coupling 
along the rungs Jr.; see Ref. IT^ Moreover, there are a 
number of further candidates for organic spin ladder ma- 
terials in the strong-coupling limit, see, e.g., Ref. El As 
an example for an inorganic system, we mention CaV205 
for which a ratio of Jl/ Jr ~ 0.1 is discussedi^. 

In the literature, bond randomness in spin ladder sys- 
tems has been studied both in the weak and strong dis- 
order limit using the real-space renormalization group 
methodiSiii, bosonizationi^, and a mapping on random- 
mass Dirac fermionsi^. Most of these studies have fo- 
cused on the stability of the ground state and the gap 
against disorder and they find that disordered spin lad- 
ders exhibit nonuniversal thermodynamic properties (see, 
e.g., Ref. ITtI) similar to disordered dimerized spin-1/2 
chainsSS. 

The plan of the paper is the following. First, we in- 
troduce the model and perform a projection on the one 
triplet subspace in Sec.^ This approximation provides 
results which are correct in leading order of the inter-rung 
coupling and are quantitatively relevant for spin ladder 
materials in the strong-coupling limit. Second, the single- 
impurity problem is solved analytically in Sec. Illll in the 
strong-coupling limit. Also, we analyze small clusters 
of bond impurities on neighboring bonds. The results, 
i.e., the eigenenergies of single- impurity (anti-) bound 
states are then compared to those of a Lanczos study 
for the full spin ladder model. In Sec. II VI finite concen- 
trations of bond-impurities are considered in the strong- 
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coupling limit. This problem is tackled both numerically 
and analytically by means of diagrammatic approaches 
[low-concentration approximation and coherent-potential 
approximation (CPA)]. The comparison of the numerical 
results with the analytical approaches provides insight 
into the validity of the latter methods. Finally, our con- 
clusions are summarized in Sec. 

II. MODEL 

The Hamiltonian of the pure two-leg spin ladder reads 

N 

Hq = ''^[JrSi^i-Si^2 + Jl0i,i-Si+i^i+Si^2-Si+i^2)]- (1) 
1=1 

^i,i{2) are spin-1/2 operators acting on site I on leg 1(2) 
and N is the number of rungs. For the remainder of the 
paper, we set Jr = 1 in all explicit computations, but we 
keep Jr in the equations for clarity. We will discuss the 
following situations 

H = Ha + H'- H' =J2 (2) 

n=l 

where A'imp is the number of modified couplings and ft,„ 
is the local perturbation at site In leading to either a 
modified on-site rung interaction = Jr -|- 5 Jr or a 
modified leg coupling J^ = Jl -|- (5Jl, connecting sites In 
and /„ + 1. Explicitly, /i„ reads 

hn = '^JrS';^,! • 5i„,2 (3) 

K = 5.hSi^^j-Si^^,^j- i = l,2. (4) 

The effect of modified interactions on the one-triplet dis- 
persion will be discussed in the strong-coupling limit 
Jl ^ Jr by projecting on the one-triplet subspace. 
Therefore, all terms contained in Hq destroying or creat- 
ing two triplet excitations are neglected. For the latter 
application of diagrammatic techniques, it is useful to 
map the spin operators on so called bond-operatorsSLSS 
sl^\t^^\] a = x,y, z. sj creates a singlet on the Ith rung 

out of the vacuum state | 0) and t]^ ^ creates a triplet 
excitation with orientation a, respectively. The exact 
representation of S'j'^2)^ — 2;,?/, z, in terms of bond- 
operators reads^^ 

S^^^ = (l/2){±4t„,, ± tl.si - le^p^tl i^^i) . (5) 

The plus sign corresponds to j = 1 and the minus sign 
to j = 2; j labeling the leg. To avoid unphysical double 
occupancies one has to impose the local constraint (sum- 
mation over repeated indices is implied in the following) 

4^l+tUta,l^^- (6) 



Projecting on the one-triplet subspace and thereby ap- 
plying a Holstein-Primakoff type of approximationSLSi 
Si ~ s\ 1 results in the effective Hamiltonian 

I I 

where we have dropped irrelevant additive constants. 
Jfo,cff is diagonalized by a Fourier transformation t^^ j = 
(l/^/]V) Y.k e^^'^'tli leading to 

Ho,e« = ^ktl.t^^k (8) 

k 

and the dispersion relation of one-triplet excitations iam 

Efc = Jr + JLCOs(fc). (9) 

The perturbations /i„ caused by modifications of the 
exchange couplings are expressed in terms of bond- 
operators as follows 

k.ki 

with the potentials fR(L)(^7^i) given by 

vnik,h) = (5JRe^'"^'= (11) 

where Afc ^ ki — k is the momentum transferred in a 
scattering process. All together, the effective Hamilto- 
nian takes the form 

-f^cff = J^o.cft + H'^fi; H^g = ^ hn- (13) 

III. THE SINGLE-IMPURITY PROBLEM AND 
SMALL IMPURITY CLUSTERS 

The solution of the one-impurity problem for the effec- 
tive Hamiltonian Eq. (|13|1 . i.e., one modified coupling, is 
derived from Schrodinger's equation in real space. Here, 
we briefly outline the procedure and our results, refer- 
ring the reader to the literatureSi for details. In addition, 
small impurity clusters are addressed and we study the 
problem of one impurity in the full model Eq. (0) using 
the Lanczos method. 



A. Single-impurity problem in real space 

Schrodinger's equation can be cast in the form 

[l-G'{E)H'^^]\i^)^0, (14) 



3 



where G^{E) = {E — Hq,cs)~^ is the free Green's func- 
tion operator associated to the effective one-particle 
Hamiltonian -ffo,off in Eq. Q and | ■0) is an eigenstate 
of the Hamiltonian iJcff • A real-space representation can 
be given using the one-triplet basis where matrix ele- 
ments of G" {E) are diagonal in a and depend only on the 
distance Al = \l — l'\. In the continuum limit N ^ oo, 
they are given by 



[G'iE)] 



a,/3 

a; 



1 

2^ 



cos(fcAZ) 



E — Jfj — Jl cos k 



dk , (15) 



with a, /? = X, y, z. Although E needs to be analytically 
continued to the complex plane in order to obtain the 
retarded Green's function by setting E E + iO+, no- 
tice that for real E^ \E — JpjJ > Jl, matrix elements in 
Eq. (|15|l are real. For simplicity, we place the modified 
rung-coupling on site I = and the modified leg-coupling 
between sites I = and ^ = 1. The (anti-) bound states 
arc found by setting the determinant of / — G'^ {E) H'^^ 
to zero, and, due to the impurity location, at most the 
upper 2x2 submatrix needs to be considered. Notice, 
however, that the eigenvalues are threefold degenerate 
because of the three triplet modes a = x,y,z. 

The eigenenergy of the (anti-) bound state in the single 
rung- impurity case 5Jr ^ 0, 5Jl = is obtained as 



El, 



Jr± 



SJn ^ 0. (16) 



The plus(minus) sign in Eq. H16(l corresponds to SJr > 
((5Jr < 0). Therefore, a bound state in the spin gap, i.e., 
below the original one-triplet band, appears for SJr < 0. 
Conversely, for SJr > 0, there is an anti-bound state 
above the one-triplet band. 

Analogously, one finds the eigenenergies of (anti-) 
bound states in the single leg-impurity case for SJr = 
0, SJl 7^ 0. For SJh > 0, there are always both a bound 
and an anti-bound state, their energies given by 



Ei^h = Jr^ 



Jl 



5Jl 
2 



1 -I- ^ 



> 0. (17) 



On the other hand, we note that there are no states out- 
side the one-triplet band for — 4Jl < (SJl < 0; instead, 
we expect the appearance of resonant modes inside the 
band [see Sec. lIVBl Fig.El(c)]. Finally, for strong ferro- 
magnetic couphng SJl < — 4Jl, Eq. (fTTjl again has two 
solutions; however, we will restrict the discussion to the 
case of antiferromagnetic couplings. 

The wave function for (anti-) bound states V'q(0 '^^^ 
also be derived in a closed form: 

SJr^O ; (5Jl = 0: 

Ml)'^[G'{Ei,R.)]r (/>0), (18) 
SJr = ; (SJl^O: 

^,(0 (X [G'{E,,^)r + [G°(i?ix)]ri 

{l>l). (19) 



The width of ?/'a(0 ^''^^l space only depends on the 
ratio of (5Jr/Jl (or (5Jl/Jl, respectively). The spatial 
extent of |V'q(OP is the narrower, the larger this ratio 
is. For instance, ]?/'„(; = 4)|V|?/'„(/ = < 0.01 for 
SJr/Jl = 1,(5Jl = 0, while |^<,(/ = 4)|VlV'a(^ = 1)1' ~ 
0.15 for SJr/Jl = 1/3. 



B. Impurity clusters 

Solving Schrodinger's equation Eq. (|14|1 in real space 
allows for the discussion of small impurity clusters. We 
consider the presence of impurities of the same type lo- 
cated on some of the first Nc ladder sites. The (anti- 
) bound eigenenergies depend in principle on both the 
number of modified couplings and their distance as well 
as the perturbation (5 Jr[l] itself. It is natural to expect 
the one-impurity eigenstates to interfere when the single 
impurities come close enough. 

As before, the eigenenergies of the cluster are evalu- 
ated by setting the determinant of / — G" (E) H'^g to 
zero. Now, only the upper Nc x Nc submatrix needs to 
be considered. As an example, we give the analytical 
expression for the case of two modified rung couplings 
on neighboring sites (Nc = 2). One solution exists for 
both SJr < (inside the gap) or SJr > (above the 
one-triplet band). Their eigenenergies E2.R read 



E' 



2,R 



Jl Jr ±Jl + 2 Jl(5Jr ± 2(5Jr( Jr. + SJj, 
Jl ± 2SJr 



(20) 



The plus sign has to be used for SJr > and the minus 
sign in the opposite case. We have computed similar 
expressions for clusters up to Nc = 5. In Sec. IIV CI 
we will show that the influence of such clusters explains 
the details of the peak structure in the density of states 
obtained by numerical diagonalization of systems with a 
finite impurity concentration. 



C. Comparison with exact diagonalization 

To test the region of validity of the results derived 
above to first order in Jl/Jr, we now compare them 
to numerical results for the full spin ladder Hamiltonian 
Eq. with one modified rung or leg coupling. 

We have exploited S'^-conservation, spin-inversion 
symmetry, reflection symmetry at the impurity bond 
and, in the case of a rung impurity, exchange symmetry 
of both legs. Usually, one uses periodic boundary con- 
ditions and exploits translational invariance. Although 
the latter is not possible if an impurity is present, we 
still apply periodic boundary conditions along the legs 
in order to minimize surface effects. Finite-size effects 
turn out to be smallest for an even number of rungs. We 
therefore concentrate on systems with = 4, 6, 8, 10, 
12, and 14 rungs. The largest dimension is slightly above 
10 million and occurs for = 14 rungs (28 spins) and 
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one leg impurity (where exchange symmetry of the legs 
is absent). 

In each of the relevant subspaccs we have com- 
puted the lowest eigenvalue using the Lanczos proce- 
dure. The results for the lowest excitation energy E at 
finite N have then been extrapolated to the thermody- 
namic limit N ^ oo using the Vanden-Broeck-Schwartz- 
algorithm2Si2Sr with a = —1. For the pure ladder this 
yields estimates for the spin gap shown by the open 
circles in Fig. Using TV < 14, we find a value of 
(0.5025 ± 0.0008) Jr at Jl = Jr, in excellent agreement 
with accepted values for this case (see section III. A of 
Ref. ^3 for a summary). As a further comparison, the 
full line in Fig. ^ shows a [7,6] Fade approximant to the 
13th order series for the spin gap of the pure ladder of 
Ref. U 

A finite number of impurities (vanishing density) does 
not affect the one-triplet band in the thermodynamic 
limit. Hence, the result for the spin gap in the pure case 
also corresponds to the lower boundary of the one-triplet 
band if impurities are present. 

Turning now to the case of one impurity, we concen- 
trate on those situations where we may expect the lowest 
excitation to be a bound state at the impurity, namely 
(5Jr < or SJh > 0, respectively. To understand the 
finite-size behavior of systems with one impurity, it is 
important to realize that there are now two competing 
length scales involved. On the one hand, there is the cor- 
relation length of the pure system, and on the other hand, 
the spatial extent of the impurity wave function needs to 
be considered [see Eqs. l(TH)l and ((T^ ]. Indeed, the typi- 
cal width of the impurity wave function, which depends 
on the actual choice of parameters, can be (much) big- 
ger than the correlation length. This interplay leads to a 
crossover in the finite-size behavior. 

For small Jl , the energy of the impurity level increases 
with system size which is in contrast to the behavior of 
the pure system where the energy of the lower edge of 
the band decreases with system size. Since the latter 
finite-size behavior is preserved at large Jl j finite-size ef- 
fects are nonmonotonic in the intermediate region, i.e., 
the finite-size behavior changes at a characteristic system 
size that increases as the impurity level approaches the 
one-triplet band. One now has to be more careful with 
the extrapolation, and we can use only those system sizes 
which arc in the asymptotic regime for large N. Accord- 
ingly, Fig. ^ shows extrapolated data points for the im- 
purity level only in a restricted region of Jl. When the 
impurity level approaches the one-triplet band, error bars 
become large, making it difficult to decide whether this 
level merges into the band or approaches it only asymp- 
totically. In any case, the numerical data demonstrate 
the presence of an impurity bound state in a wide pa- 
rameter region. 

In the limit of small Jl Jr, we can compare 
to Eqs. (|16|) and (|17|) . respectively. One indeed ob- 
serves quantitative agreement for sufficient small Jl, see 
Fig.n(a) and the inset of Fig.n](b). Note that due to the 
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SJ^ = -0.1 
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8J^ = -0.2 
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SJ^ = -0.3 




T 


5Jj, = -0.4 






Eq. (16) 






FIG. 1: Excitation energy E of the lowest level for the full 
spin ladder Hamiltonian (|5J. Panel (a): one rung impurity 
Jr = Jr + SJr. Panel (b): one leg impurity Jl = Jl -I- (5Jl. 
In all cases, the normalization is fixed to Jr = 1. Symbols 
are obtained by extrapolation of Lanczos diagonalization on 
finite systems. Open circles are for the pure system {SJr = 
and (5Jl = 0) and correspond to the spin gap; the solid 
line is a [7,6] Pade approximant to the 13th order strong- 
coupling series- for the spin gap of the pure ladder. Dashed 
lines display the analytical result for the position of the bound 
state in the effective Hamiltonian Eq. 1131 . namely Eq. (11611 
[panel (a)] and Eq. llTTl [panel (b)]. 



normalization of Fig. ^ i.e., (5Jl ^ Jl, Eq. H17|l results in 
straight lines which start a,t E = 1 for Jl = 0. At larger 
Jl deviations can be observed in Fig. ^ but an impu- 
rity level can still be seen. Accordingly, the first-order 
approximation can still be expected to be qualitatively 
correct even in a parameter region where it is no longer 
quantitatively accurate. Hence, we may use the first- 
order approximation to study several impurities and even 
finite densities which is no longer systematically possible 
by Lanczos diagonalization of the full ladder Hamilto- 
nian. 
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FIG. 2: Sketch of the diagrammatic expansion of the self 
energy E(_E) in the low-concentration hmit. G1{E) = 1/(E — 
ek) is the free one-triplet Green's function. 



IV. FINITE CONCENTRATIONS 

In this section, we discuss the effect of a finite concen- 
tration of modified couplings Jr or Jl on the one-triplet 
dispersion in the strong-coupling limit Jl ^ Jr. To 
treat this problem, we apply diagrammatic techniques, 
namely, a low-concentration approximation (LCA) and 
the coherent-potential approximation (CPA), and numer- 
ical diagonalization of large systems. We also use our 
analytical results for the eigenenergies of small impurity 
clusters to explain the details in the numerical results for 
the density of states, both for (anti-) bound states and 
resonance modes. 

Before we turn to the discussion of the analytical meth- 
ods and compare the results to those from numerical im- 
purity averaging (NAV), let us consider certain limiting 
cases. In the following, c denotes the concentration of 
impurities. Note that in the case of impurities on the 
legs, we set c = 1 if all 2N couplings are modified. 

The limiting cases are: (i) the pure system (c = 0); 
(ii) the single- impurity case (see Sec. IIII A(l : (iii) the 
case c = 1 where all couplings are equal to Jr + (5Jr 
or Jl -|- i5Jl, respectively. In the latter case and for 
i5Jr ^ 0, (5Jl = 0, a one-triplet band with the disper- 
sion Ek ~ {Jr + 6Jn) + JhCOs{k) will result, i.e., its 
center is shifted by 6Jr with respect to the center of the 
original band for c — 0. Therefore, the single-impurity 
(anti-) bound state should develop into a dispersive band 
as the concentration increases while the center e(c) of the 
band lies between Jr — \/ Jl + (<^Jr)^ < e(c) < Jr -|-(5Jr 
for SJr < 0. An analogous scenario arises for JJr > 0. 



For (5 Jr = 0, (5 Jl > 0, the triplet dispersion in the limit 
of c = 1, i.e., all Jl modified, reads Ek = Jr. + (Jl + 
(5JL)cos(fc). Thus, the bound- and anti-bound states 
appear symmetrically with respect to the center of the 
original band. On increasing the concentration c, ad- 
ditional impurity levels will appear and eventually, they 
will merge in the original band. Finally, there will be one 
broadened band possessing a bandwidth of (Jl + (5Jl)- 



A. Low-concentration and coherent-potential 
approximation 

Based on a diagrammatic expansion of the one-triplet 
Green's function in the presence of impurities, a number 



of useful methods exist to get approximate results for 
the self-energy S(£'). First, we briefly comment on the 
low-concentration approximation and second, we discuss 
results from the coherent-potential approximation. 

As in Sec. Illll we will concentrate on the single-triplet 
subspace. Thus, apart from integrating out the singlet, 
the hard-core constraint Eq. @ is automatically satisfied 
within our approximation, i.e., first-order perturbation 
theory in Jl/ Jr. 

Low- concentration approximation - Using standard 
impurity-averaging techniques (see, e.g., Ref. l2^ . the 
self-energy of the one-triplet Green's function can be 
obtained in first order in the impurity concentration. 
As the averaging procedure restores translational in- 
variance, the one-triplet Green's function Gk{E) can be 
written in terms of the Dyson equation 



Gk{E) 



1 



E-ek- ci:{E) 



(21) 



Keeping only terms linear in c implies that the self- 
energy S(i5) is equal to the T-matrix of the one- impurity 
problem. The diagrammatic expansion of the self energy 
is sketched in Fig. |21 Note that all quantities in Eq. H21|) 
become 2x2 matrices if the modified coupling connects 
two sites as realized by a leg impurity. For a more 
detailed discussion of this technique, the reader is 
referred to, e.g., Ref. [2^ 

The spectral function Ak{E) = -{l/-n)luiGk{E) is 
plotted in Fig.Elfor (a) Jl = 0.1, (5Jr = -0.1, c = 0.01 
and (b) for Jl = 0.1, (5Jl = 0.5, c = 0.01. In accordance 
with our previous results we find one bound state in 
case (a) and a bound and an anti-bound state in case 
(b). Figure 13 further reveals that first, the impurity 
levels have developed a small dispersion and second, the 
spectral weight is concentrated around k = tt for the 
bound states while it vanishes in the center of the zone, 
and vice- versa for the anti-bound states. 

Coherent-potential approximation - The coherent- 
potential approximation allows one to interpolate 
between the two limits of c = and c = 1. Here, we 
apply this method to the case of (5 Jr < 0, (5Jl = 0. The 
self-energy is obtained from a self-consistent solution of 
the equation^ 



■SiE) 



c6Jn 



1 - G{E) [<5Jr - Y.{E)] ■ 



(22) 



Rather than deriving this equation (see Ref. |23 for de- 
tails), let us mention some features of this method: (i) 
the self-energy is symmetric under exchange of host and 
impurity sites, i.e., c and 1 — c and the respective replace- 
ment of the coupling constants; and (ii) it gives qualita- 
tively correct results for the density of states for inter- 
mediate concentrations. We note that in contrast to the 
low-concentration approximation [see, e.g.. Fig. ^ (a)], 
the CPA does not lead to a sharp peak in the density 
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FIG. 3: Spectral function Ak{E) in low-concentration ap- 
proximation (LCA) for (a) perturbed rung couplings with 
(5Jr = —0.1, c = 0.01; (b) perturbed leg couplings with 
(5Jl = 0.5, c = 0.01. The dashed lines mark the positions 
of the (anti-) bound states from Eqs. (1161 and 1171 1. Jr = 1 
and Jl =0.1 in both cases. 



of states at the position of the impurity level even for 
low concentrations. This can, for example, be seen in 
Fig- E (b) for c = 0.1. We have, however, checked that 
in both diagrammatic approaches, the total weight in the 
impurity levels is the same and that it grows linearly with 
the impurity concentration, as expected. 



B. Numerical results and comparison 

Now we compare the analytical results with a numeri- 
cal diagonalization of the effective Hamiltonian on large 
systems and sampling over several realizations at fixed 
concentration. The effective Hamiltonian Eq. 
has been diagonalized on finite systems with N = 10'^ 
rungs for different choices of impurity concentrations 
c for both types of bond impurities. The density of 
states (DOS) is obtained from binning the eigenvalues, 
the bin- width of typically Ai? ^ 10 ^''Jr determining 
the resolution in Figs. 01 to Results are shown for 
Jl = 0.1, (5Jr = -0.1 in Fig. |i|(panel (a): c = 0.01; 
(b): c = 0.1; (c): c = 0.3). Note that, according 
to Eq. H16() . the position of the single-impurity level is 
Ei,K = 0.8586Jr,. The following features are observed: 
(i) for increasing concentration, additional peaks appear 
in the vicinity of the one-impurity level. They stem from 
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Jl=0.I,8J,=-0.1,N=10" 
(a) c=0.01 

LCA 

NAV 



(b)c=0.1 
CPA 



(c) c=0.3 
CPA 




FIG. 4: Density of states (DOS) at finite concentration of 
modified couplings Jr: numerical data (NAV, solid line) 
for spin ladders with A'' = lO'^ rungs and Jr = 1, Jl = 
0.1, (5Jr = —0.1 (concentration: panel (a) c = 0.01; panel 
(b) c = 0.1; panel (c) c = 0.3.). Dashed line in panel (a): 
low-concentration approximation (LCA); in panel (b) and (c): 
CPA. 



impurity clusters, i.e., impurities occupying neighboring 
sites, as will be discussed in more detail below, (ii) The 
bound state level develops into a band centered around 
Jr -|- i5Jr = 0.9 as a function of concentration c. No- 
tice that larger concentrations c > 0.5 are conveniently 
realized by setting c — > 1 — c, Jr Jr + ^Jr,, and 
(5Jr — > — i5Jr. (iii) Inside the original band, the curve 
is not smooth, but displays small oscillations. These fea- 
tures are neither due to finite-size effects nor due to low 
statistics (the density of states has been obtained by aver- 
aging over typically a few thousand random realizations 
at fixed concentration). As we shall discuss below, their 
origin can also be related to the effect of impurity clus- 
ters. 

Let us now comment on the comparison of the numeri- 
cal with the analytical results. By integrating Ak{E) over 
the momentum k, the density of states n{E) is obtained. 
Results from the LCA are compared to the numerical im- 
purity averaging in the case of Jl = 0.1, (5Jr = —0.1, c = 
0.01 in Fig. 21 (a)- Both approaches agree well with re- 
gard to the position of the main impurity level. The 
comparison with the results from the CPA for c = 0.1 
[Fig. H (b)] and c = 0.3 [Fig. |3| (c)] shows that this ap- 
proach gives qualitatively reasonable results even at fairly 
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FIG. 5: Density of states (DOS) at finite concentration of 
modified couplings J^: numerical data (NAV, solid line) for 
spin ladders with A'^ = 10^ rungs and Jr — 1, Jl = 0.1, 
and (5Jl = 0.5 (concentration: panel (a) c = 0.01; panel 

(b) c = 0.1; panel (c) c = O.9.). Dashed line in panel (a); 
low-concentration approximation (LCA). Inset of panel (a): 
structure of the DOS in the vicinity of the bound state. Panel 

(c) : one observes resonance modes inside the band for c — 0.9; 



the case shown here is equivalent to Jl 
and c — 0.1. 



: 0.6, (5Jl 



-0.5 < 0, 



large concentrations. At c = 0.3, the impurity levels and 
the original band start to merge. 

For a finite concentration of leg couplings, the numer- 
ical results confirm our qualitative expectations. The 
data are shown in Fig. for Jl = 0.1,i5Jl = 0.5 (panel 
(a): c = 0.01; (b): c = 0.1; (c): c = 0.9). For clarity, 
we note that the possible impurity configurations are: (i) 
one modified coupling on one leg, connecting, e.g., rung 
I and ^ -|- 1; and (ii) both couplings between rung I and 
I + 1 modified. Both cases arc taken into account in the 
numerical implementation. 

The impurity levels occur symmetrically with respect 
to the center of the band. On increasing the concentra- 
tion c, the original band widens and eventually includes 
all impurity levels (see Fig. |S1 (c) for c = 0.9). The in- 
fluence of impurities is now visible as resonance modes 
inside the band. Note that the last case is equivalent to 
Jl = 0.6, 6Jl = -0.5, and c = 0.1. 

Comparing to the LCA from Eq. H21|) for the case of 
Jl = 0.1, (5 Jl = 0.5, c = 0.01 [see Fig. El (a)], we see that 
the positions of the highest peaks seen in the numerical 
data and the analytical result almost coincide, similar 
to the case shown in Fig. 0] (a). The inset of Fig. [S] (a) 
contains a zoom of the region around the lower single-leg 
impurity peak with -Bi^l = 0.811 Jr revealing the pres- 
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FIG. 6: Comparison of eigenenergies of bound states induced 
by rung- impurity clusters in a clean system, derived analyti- 
cally, with numerical diagonalization of Hcb on large systems 
with a finite concentration of modified rung couplings. The 
density of states (DOS) in the vicinity of the one-impurity 
level is shown. The patterns in parenthesis denote different 
types of impurity clusters: with 1, we indicate the relative po- 
sition of the rung impurities in the cluster sequence. The nu- 
merical data (NAV, solid line) correspond to Jr = 1, Jl = 0.1, 
(5Jr = —0.1 and c — 0.1. The letters 'a' to 'j' relate the peaks 
in the DOS to certain impurity clusters, which are listed in 
the legend. 



ence of several less pronounced structures. 

In summary, the analytical approaches give fair re- 
sults for the overall structure of the density of states 
even at large concentrations as exemplified in the case 
of 5 Jr < 0, (5Jl = 0. The effects of impurity clusters are 
not taken into account in the diagrammatic description. 



C. Analytical results for small impurity clusters 

Next we analyze the details of the peak structure of 
the impurity levels arising from a variety of small impu- 
rity clusters for the case of modified rung couplings. The 
clusters analyzed include, e.g., the patterns (11), (101), 
(111), (1001), (1101), (1111), and (10001) where, in this 
notation, 1 indicates a modified rung in the cluster se- 
quence and indicates no impurity placed on a rung site. 
For example, (11) denotes two rung impurities on neigh- 
boring sites in an otherwise clean system. 

We find that the energy eigenvalues corresponding to 
the main peaks outside the triplet band observed in the 
numerical results for the density of states (see the pre- 
vious section) can be associated with the contribution 
of certain impurity clusters. In particular, in Fig. |S1 we 
show the matching between the numerical bound state 
structure and the analytically computed eigenenergies 
corresponding to different small clusters placed in an oth- 
erwise clean system for the case of Jl = 0.1, (5Jr, = —0.1 
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FIG. 7: Matching between the peak structure inside the 
triplet band obtained numerically at finite concentration of 
rung impurities (NAV, solid line) and analytically for sys- 
tems with different types of impurity clusters (dashed line), 
for Jr = 1, Jl = 0.01, SJk = -0.05 and c = 0.1. 

and c = 0.1 (see Fig. 0| (b) for the full density of states 
in the same case). 

Let us now discuss the features inside the triplet band 
in the presence of impurities. The density of states (DOS) 
n(E) can be evaluated from 

n{E) = --ImTTGiE) = n^{E) + n,,,,^{E) , (23) 

TT 

where the free density of states^S is no{E) = 
— (I/tt) ImTrG'°(£') and the contribution from the im- 
purities riirapiE) can be written as 

nimAE) = In [Oct (/ - G° (E) iJ^ff)]. (24) 

We use Eq. H24|) to compute the contributions from 
particular impurity clusters to the DOS inside the triplet 
band. Moreover, from the amplitudes of the different 
peaks outside the one-triplet band, one can read off, at 
least qualitatively, the distribution of probabilities for 
the presence of the different clusters in a random sample 
of impurities for a given concentration. This information 
is in turn used to weight the influence of each cluster on 
the peak structure inside the triplet band. In Fig. [7| we 
show an example where the main peaks are associated 
with the corresponding cluster contributions. We are 
able to match, in this particular case, the central peak, 
which is slightly shifted to the right of the band-center. 



with the patterns (1101) and (101), and the two ones, 
which are almost symmetrical with respect to the band 
center, with a contribution from the (1001) cluster. 

This analysis explains on the one hand the appearance 
of the various localized modes and on the other hand 
it gives evidence that the peaks inside the triplet band 
originate from the existence of impurity clusters. 

V. DISCUSSION AND CONCLUSIONS 



In this paper we have studied the appearance of bound 
states in the spin gap of spin ladders with bond im- 
purities. Both the cases of modified rung and leg cou- 
plings have been considered. We have derived analytical 
results for the position of bound states in the strong- 
coupling limit equivalent to first order perturbation the- 
ory in Jl/Jr. The existence of impurity induced bound 
states has been verified by a Lanczos study of finite spin 
ladders with one impurity and < Jl < Jr. and we 
find that our analytical results are quantitatively correct 
foi' Jt^ ^ Jr/IO and that a qualitative agreement is still 
found for larger Jl ■ Recently discovered spin ladder ma- 
terials such as, for example, (C5Hi2N)2CuBr4 (Ref. IT^ 
or CaV205 (Ref. ITsh fall in this range of parameters. 

Further, we have discussed the density of states in the 
presence of a finite concentration of impurities in the 
limit of Jl ^ Jr both numerically and analytically. The 
comparison of the different approaches shows that dia- 
grammatic methods give quantitatively correct results for 
small impurity concentrations and, furthermore, a qual- 
itatively correct picture is obtained for large impurity 
concentrations. As the diagrammatic approaches neglect 
the interference of impurities and the effect of impurity 
clusters, we have presented a careful analysis of systems 
with small impurity clusters which allows us to under- 
stand details visible in the density of states. Natural 
extensions of this work, i.e., the computation of observ- 
ables and the discussion of systems with arbitrary ratios 
of Jl/Jr,, are left for future work. Nevertheless, our re- 
sults already imply the appearance of additional features 
in the spin gap which could be observed by, e.g., optical 
experiments on bond-disordered spin ladder materials. 
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